Non-linear amplification of small spin precession using long range dipolar interactions 
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In measurements of small signals using spin precession the precession angle usually grows linearly 
in time. We show that a dynamic instability caused by spin interactions can lead to an exponentially 
growing spin precession angle, amplifying small signals and raising them above the noise level of a 
detection system. We demonstrate amplification by a factor of greater than 8 of a spin precession 
signal due to a small magnetic field gradient in a spherical cell filled with hyperpolarized liquid 
129 Xe. This technique can improve the sensitivity in many measurements that are limited by the 
noise of the detection system, rather then the fundamental spin-projection noise. 
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Observation of spin precession signals forms the basis 
of such prevalent experimental techniques as NMR and 
EPR. It is also used in searches for physics beyond the 
Standard Model [J 0, IS B an d sensitive magnetome- 
tery 0. Hence, there is significant interest in the de- 
velopment of general techniques for increasing the sensi- 
tivity of spin precession measurements. Several methods 
for reducing spin-projection noise using quantum non- 
demolition measurements have been explored jy, [jj and 
it has been shown that in some cases they can lead to im- 
provements in sensitivity 0, Ej . In this Letter we demon- 
strate a different technique that increases the sensitivity 
by amplifying the spin precession signal rather than re- 
ducing the noise. 

The amplification technique is based on the exponen- 
tial growth of the spin precession angle in systems with 
a dynamic instability caused by collective spin interac- 
tions. Such instabilities can be caused by a variety of 
interactions, for example, magnetic dipolar fields in a 
nuclear-spin-polarized liquid [T(J, HH H2 or electron-spin 
polarized gas|l3|, spin-exchange collisions in an alkali- 
metal vapor |l4j or mixtures of alkali-metal and noble-gas 
atoms pjj. This amplification technique can be used in 
a search for a permanent electric dipole moment in liq- 
uid 129 Xe [3]. It is also likely to find applications in a 
variety of other systems with strong dipolar interactions, 
such as cold atomic gases 17] and polar molecules [l8j . 

Consider first an ensemble of non-interacting spins 
with a gyromagnetic ratio 7 initially polarized in the x di- 
rection and precessing in a small magnetic field B z . The 
spin precession signal (S y ) = r )(S x )B z t grows linearly in 
time for jB z t <C 1. The measurement time t m is usually 
limited by spin relaxation processes and determines, to- 
gether with the precision of spin measurements S((S y )), 
the sensitivity to the magnetic field B z 
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or any other interaction coupling to the spins. In the 
presence of a dynamic instability, the initial spin preces- 
sion away from a point of unstable equilibrium can be 
generally written as (S y ) = r y(S x )B z sinh(/3t)//3, where 
/3 is a growth rate characterizing the strength of spin in- 
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Hence, for the same uncertainty in the measurement of 
(S y ), the sensitivity to B z is improved by a factor of 
G = sinh(/?t m )//?t m . It will be shown that quantum (as 
well as non-quantum) fluctuations of (S y ) are also ampli- 
fied, so this technique cannot be used to increase the sen- 
sitivity in measurements limited by the spin-projection 
noise. However, the majority of experiments are not 
limited by quantum fluctuations. For a small number 
of spins the detector sensitivity is usually insufficient to 
measure the spin-projection noise of TV 1 / 2 spins, while 
for a large number of particles the dynamic range of the 
measurement system is often insufficient to measure a 
signal with a fractional uncertainty of A/" -1 / 2 . Ampli- 
fying the spin-precession signal before detection reduces 
the requirements for both the sensitivity and the dynamic 
range of the measurement system. Optical methods al- 
low efficient detection of electron spins and some nuclear 
spins 3] in atoms or molecules with convenient optical 
transitions. However, for the majority of nuclei optical 
detection methods are not practical and magnetic detec- 
tion, using, for example, magnetic resonance force mi- 
croscopy, has not yet reached the sensitivity where it is 
limited by the spin projection noise [HI [2(| . Therefore, 
non-linear amplification can lead to particularly large im- 
provements in precision measurements relying on nuclear 
spin precession. 

Here we use long-range magnetic dipolar interactions 
between nuclear spins that lead to exponential amplifi- 
cation of sp in precession due to a magnetic field gradient 
[ill Ufl l21j . It has also been shown that long-range dipo- 
lar fields in conjunction with radiation damping due to 
coupling with an NMR coil lead to an increased sensitiv- 
ity to initial conditions and chaos [22( • To amplify a small 
spin precession signal above detector noise it is important 
that the dynamic instability involves only spin interac- 
tions, since instabilities caused by the feedback from the 
detection system would couple the detector noise, such 
as the Johnson noise of the NMR coil, back to the spins. 
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FIG. 1: Low field NMR setup (view from above). Polarized 
liquid 129 Xe is contained in a spherical cell maintained at 
173K by flowing N2 gas through a vacuum insulated column. 
High-T c SQUIDs are submerged in LN2 contained in a glass 
dewar. Inset: configuration of the SQUIDs, applied magnetic 
field, the magnetization, and the rotatable membrane. 



We measure spin precession using SQUID magnetome- 
ters that do not have a significant back-reaction on the 
spins and show that under well controlled experimental 
conditions the dynamic instability due to collective spin 
interactions can be used to amplify small spin precession 
signals in a predictable way. 

Our measurements are performed in a spherical cell 
containing hyperpolarized liquid 129 Xe (Fig. 1). Liquid 
129 Xe has a remarkably long spin relaxation time and 
the spin dynamics is dominated by the effects of long- 
range magnetic dipolar fields. In the spherical geometry 
an analytic solution can be found using a perturbation 
expansion in a nearly uniform magnetic field Ho 
We are primarily interested in the first-order longitudinal 
magnetic field gradient g, H = (Ho + gz)z, but will also 
consider other magnetic field gradients which inevitably 
arise due to experimental imperfections. For longitudinal 
gradients that preserve cylindrical symmetry the magne- 
tization profile can be expanded in a Taylor series, 
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where R is the radius of the cell. Only gradients of the 
magnetization create dipolar fields in a spherical cell, for 
example, a linear magnetization gradient m^ 1 ' ) creates 
only a linear dipolar magnetic field, which, in the rotating 
frame, is given by 
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The time evolution of the magnetization is determined 
by the Bloch equations dM/dt = 7M x (B d + H). If the 
magnetization is nearly uniform, m^' fe ) <C 1, they can 
be reduced to a system of linear first-order differential 
equations for m^ hk \ 

We consider first the simplest case when only the linear 
field gradient g is present and the initial uniform magne- 
tization Mo is tipped into the x direction of the rotating 



frame by a tt/2 pulse. Substituting Eqns. (|3J) and (O into 
the Bloch equations we find that only linear magnetiza- 
tion gradients grow as long as m^'^ <C 1, in particular, 
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Here f3 is proportional to the strength of the long-range 
dipolar interactions. We measure m^ 1 ' ^ experimentally 
by placing two SQUID detectors near the spherical cell 
as illustrated in Fig. Q an d measuring the phase differ- 
ence Acj) between the NMR signals induced in the two 
SQUIDs. For small m^°\ A</> = (m^ ,0 \ where £ is a 
numerical factor that depends on the geometry, for our 
dimensions £ = 0.46 ± 0.01. Thus, the phase difference 
A(/> is proportional to the applied magnetic field gradient 
g and grows exponentially in time, increasing the sensi- 
tivity to g by a factor G = smh((3i) / (3t. For M = 100 
/iG, which is easy to realize experimentally with hyper- 
polarized 129 Xe, (3 = 1.75 sec -1 , so that a very large 
amplification factor can be achieved in a short time, for 
example G = 360 after 5 seconds. 

One of the main challenges to realizing such high gains 
is to achieve sufficient control over the initial conditions 
and non-linear evolution of the system, so that the dy- 
namic instability gives rise to a phase difference Ac/) that 
remains proportional to g even in the presence of various 
experimental imperfections. We developed a set of nu- 
merical and analytical methods for analyzing these effects 
[2^ . Since our goal is to achieve very high sensitivity to 
a small first-order longitudinal magnetic field gradient g, 
we generally assume that it is smaller than other gradi- 
ents that are not measured directly. We find that the 
presence of transverse or higher order longitudinal gra- 
dients as well as initial magnetization inhomogeneities 
cause an abrupt non-linear decay of the overall magne- 
tization. The time until the decay t c depends on the 
size of the inhomogeneities relative to Mo and limits the 
achievable gain to smh((3t c ) / (3t c . Inhomogeneities of the 
applied field symmetric with respect to the z direction 
do not change the evolution of A</>, which remains pro- 
portional to g until the collapse of the magnetization, 
as shown in Fig. [2^ Higher order z-odd longitudinal 
gradients do generate a phase difference (Fig. 03). How- 
ever, the contributions of different magnetic field gra- 
dients to the phase difference add linearly as long as 
m (*>*0 1 and the effects of higher order odd gradients 
can be subtracted if they remain constant, as illustrated 
in Fig. I2t>- While higher order magnetization gradients 
can grow with a time constant up to 2.5 times faster than 
the first-order gradient, it can be shown using a pertur- 
bation expansion that the first moment of the magneti- 
zation d = J zMydV always grows with an exponential 
constant given by Eq. and is proportional to the first 
moment of the magnetic field b = J zB z dV . The phase 
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FIG. 2: Numerical simulations |23| of the SQUID signal (left 
axis) and the phase difference between SQUIDs (right axis) 
for Mo = 100 /xG and a small longitudinal field gradient 
g = 0.1/xG/cm (solid lines), a) An additional larger trans- 
verse gradient g± = 2 /xG/cm (dashed line) or a second-order 
longitudinal gradient gi — 1 /xG/cm 2 (dash-dot) do not affect 
the phase difference until the SQUID signal begins to decay, 
b) Effects of an additional z-odd third-order longitudinal gra- 
dient #3 = 0.8/xG/cm 3 (squares). Stars show the phase evo- 
lution in the presence of gs but for g — 0. The difference 
between the phase for g = 0.1 /xG/cm and g = (triangles) 
follows the solid line corresponding to the pure linear gradient 
g until the magnetization begins to collapse. The third-order 
gradient generates a background phase that can be subtracted 
to determine a change in g between successive measurements. 



difference between the SQUID signals is approximately 
proportional to the first moment of the magnetization d 
and is not significantly affected by the growth of higher 
order gradients. For example, in Fig.Eb) the overall sig- 
nal decays at about 3 sec due to large first and third-order 
magnetization gradients but the phase difference re- 
mains much less than 1. 

Hence, the phase difference A(j) can be used to mea- 
sure a very small linear gradient g in the presence of 
larger inhomogeneities if all magnetic field and magneti- 
zation inhomogeneities are much smaller than Mq. The 
ultimate sensitivity is limited by fluctuations of the gra- 
dients between successive measurements. In addition to 
fluctuations of g, which is the quantity being measured, 
the phase difference will be affected by the fluctuations 
in the initial magnetization gradients m^ 1 ' ^ and m^ 1 ' ^ 
and, to a smaller degree, higher order z-odd gradients 
of the magnetic field and the magnetization. In par- 
ticular, fluctuations of m^ 1 ' ^ and m^ 1 ' ^, either due to 
spin-projection noise or experimental imperfections, set a 
limit on the magnetic field gradient sensitivity on the or- 
der of Sg = &TT\/2M()5m£^ /IbR and similar for 5rriz^ . 
The shot noise fluctuations of 129 Xe magnetization gener- 
ate a magnetic field gradient on the order of 10 -13 G/cm. 

Hyperpolarized 129 Xe is produced using the standard 
method of spin exchange optical pumping [3, El • The 
polarized gas is condensed in a spherical glass cell held at 
173 K as shown in Fig.^ The cell, with an inner radius 
R = 0.55 cm, is constructed from two concave hemispher- 
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FIG. 3: a) Oscillating transverse magnetization following a 
7r/2 pulse. After the signal drops to 90% of its initial value 
a second pulse is applied to realign the magnetization with 
the longitudinal direction, b) Phase difference between the 
SQUID signals. Overlaying the data (dashed line) is a fit 
based on Eq. (0. The dash-dot line is the expected phase 
evolution in the non-interacting case, illustrating that the sig- 
nal would barely be detectable. 



ical lenses glued together with UV curing cement. Inside 
the cell is an octagonal silicon membrane 25 /xm thick, 
with a diameter of 1.05 cm. The membrane is connected 
to a stepper motor outside the magnetic shields via a 0.2 
mm glass wire to mix the sample, ensuring uniformity of 
the polarization. In addition to mixing the sample, the 
membrane inhibits convection across the cell due to small 
temperature gradients which can wash out the longitu- 
dinal gradient of the magnetization. A set of coils inside 
the shields create a 10 mG uniform magnetic field and 
allow application of RF pulses and control of linear and 
quadratic magnetic field gradients. The NMR signal is 
detected using high-T c SQUID detectors. The pick-up 
coil of each SQUID detector is an 8x8 mm square loop 
located approximately 1.6 cm from the center of the cell 
and tilted by ±45° relative to the magnetic field. 

In our experimental system, the time scale of the dipo- 
lar interactions is much smaller than the spin relaxation 
time or the time needed to polarize a fresh sample of 
129 Xe. In order to make multiple measurements on a sin- 
gle sample of polarized xenon, we first apply a tt/2 pulse 
and monitor in real time the SQUID signals. When the 
NMR signal drops to 90% of its initial value, a second 
7r/2 pulse is applied, realigning the magnetization with 
the holding field. The silicon membrane is then oscil- 
lated back and forth to erase the magnetization inhomo- 
geneities developed in the previous trial. 

Fig.[2K) shows the oscillating transverse magnetization 
and Fig. |3t>) shows the phase difference between the two 
SQUID signals. We determine the value of (3 from the 
magnitude of the NMR signal and fit the phase differ- 
ence to Eq. © with g as the only free parameter. The 
dash-dot line shows the expected evolution of the phase 
difference for the same gradient in the absence of dipo- 
lar interactions, demonstrating that without amplifica- 
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FIG. 4: a) Measurement of a small gradient g alternated be- 
tween successive trials. Stars show the applied linear gradient, 
squares show the gradient measured using non-linear spin evo- 
lution, b) Gain G associated with non-linear spin evolution. 
The gain drops when the sample is not mixed in the shaded 
region, demonstrating the significance of initial magnetization 
inhomogeneit ies . 



tion the phase difference would be barely above the noise 
level of the detection system. For this measurement the 
phase is amplified by a factor of 9.5 before the magneti- 
zation drops to 90% of its initial value. 

By applying a series of double tt/2 pulses we can make 
repeated measurements of the magnetic field gradient. 
Fig. SJl) shows data where the applied longitudinal gra- 
dient is oscillated with an amplitude of 1 /xG/cm between 
trials. The stars show the applied gradient, the squares 
show the gradient measured by the non-linear spin evolu- 
tion, indicating that the amplified signal follows the ap- 



plied gradient. Slight differences between the two curves 
are due to noise in the magnetic field gradient as well 
as possible imperfections in the erasing of magnetization 
gradients between successive trials. Fig. EJd) shows the 
gain parameter for the same data set. We associate the 
rising gain at the beginning of the data set with a decay of 
the magnetization inhomogeneities developed during col- 
lection of liquid 129 Xe in the cell. In the shaded region 
of the plot we did not mix the magnetization with the 
membrane before the measurement, resulting in a drop 
of the gain as well. Numerical simulations indicate that 
the gain is likely limited by higher order gradients, for 
example a second-order magnetic field gradient on the 
order of 1 /iG/cm 2 , which can not be excluded based on 
our mapping of ambient fields, is sufficient to limit the 
gain to about 10. 

In conclusion, we have demonstrated that non-linear 
dynamics arising from long range dipolar interactions can 
be used to amplify small spin precession signals, improv- 
ing the signal-to-noise ratio under conditions where lim- 
itations of the spin detection system dominate the spin 
projection noise. By amplifying the signal before detec- 
tion, this technique reduces the requirements on the sen- 
sitivity of the detection technique as well as its dynamic 
range. In addition to precision measurements, this tech- 
nique can potentially be used to amplify small spin pre- 
cession signals in various MRI applications, allowing, for 
example, direct detection and imaging of the magnetic 
fields generated by neurons with MRI |25j |. Initial inho- 
mogeneities of the magnetization, caused, for example, 
by very slight differences of T\ in tissues, can also be am- 
plified. We thank DOE, NSF, the Packard Foundation 
and Princeton University for support of this project. 
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